Yang- Yang thermometry and momentum distribution of a trapped one-dimensional Bose gas 
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We describe the use of the exact Yang- Yang solutions for the one-dimensional Bose gas to enable accurate 
kinetic-energy thermometry based on the root-mean-square width of an experimentally measured momentum 
distribution. Furthermore, we use the stochastic projected Gross-Pitaevskii theory to provide the first quanti- 
tative description of the full momentum distribution measurements of Van Amerongen et al., Phys. Rev. Lett. 
100, 090402 (2008). We find the fitted temperatures from the stochastic projected Gross-Pitaevskii approach 
are in excellent agreement with those determined by Yang- Yang kinetic-energy thermometry. 
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Ultracold gases offer a unique opportunity to study fun- 
damental problems in quantum many-body physics, allowing 
experimental observations to be compared directly with mi- 
croscopic theories. An area of significant recent interest has 
been the measurement of thermodynamic relations . The 
one-dimensional (ID) Bose gas with repulsive interactions 
has emerged as a paradigm system because exact solutions are 
available for both eigenstates J3l and thermodynamic quanti- 
ties (see Ref. |0] and references therein). Furthermore, this 
system exhibits a surprisingly rich variety of regimes J3, 01 
connected by broad crossovers. Most studies of the ID Bose 
gas have focused on the position-space distributions HIM] 
and local correlations J3[Il[IIl-[9> which can be directly ob- 
tained from the exact theories. 

A recent experiment by Van Amerongen et al. ifioll mea- 
sured the position and momentum distributions of a trapped 
ID Bose gas throughout the crossover from an ideal gas to the 
quasicondensate regime. The position-space measurements 
were compared with the Yang-Yang (YY) thermodynamic so- 
lutions [5] within the local density approximation (LDA), and 
showed smooth behavior throughout the crossover. In con- 
trast, the momentum distributions showed a pronounced tem- 
perature dependence, and have been unexplained by theory to 
date. Previous work on the momentu m proper ties of the ID 



Bose gas has focused on limiting cases IU2U13U15I1 



Here we investigate the momentum properties of the ID 
Bose gas and their application to thermometry through mea- 
surements of the system kinetic energy. Our methods provide 
a reliable foundation for accurate thermometry in all regimes 
of a ID Bose gas with repulsive interactions, including the 
strongly-correlated regime. This approach is reminiscent of 
molecular dynamics calculations, where the average kinetic 
energy per particle is a direct measure of the temperature lfl6ll . 

First, we use the exact YY thermodynamic formalism J^t] 
to calculate the root-mean-square (rms) width of the momen- 
tum distribution, which is equivalent to determining the av- 
erage kinetic energy per particle. In combination with the 
LDA for trapped (nonuniform) quasi- ID systems, we show 
how the YY kinetic energy results can be applied to accu- 



rate thermometry for a broad range of conditions that are rel- 
evant to ongoing experimental and theoretical efforts. We 
refer to this approach as YY thermometry. Second, we 
present the first quantitative calculation of the full momen- 
tum distribution for a trapped quasi- ID Bose gas using the 
stochastic projected Gross-Pitaevskii equation (SPGPE) tech- 
nique Il7ll . (Calculations of the position-space distribution 
based on a related formalism IU8I1 were recently described 
in lfl9ll .) We find excellent agreement between the SPGPE 
results and the momentum-space measurements reported in 
Ref. d. Finally, we compare the kinetic energy and tem- 
perature predictions of the YY thermometry with those of the 
SPGPE and also find excellent agreement between the two 
approaches. While the SPGPE technique is limited to the de- 
generate yet high-temperature weakly interacting regime, the 
YY thermodynamic formalism applies to all repulsive inter- 
action strengths and temperatures and therefore its usefulness 
extends beyond the regimes studied here. 

A uniform ID Bose gas in the thermodynamic limit is com- 
pletely characterized by two parameters Jj-ISI] ; the dimension- 
less interaction strength 7 = mg/h 2 p and the reduced tem- 
perature t = 2kBTH 2 /mg 2 , where p is the linear (ID) den- 
sity, g ~ 2fku±a is the effective ID coupling strength lEoll . 
a is the 3D s-wave scattering length, and uj± the transverse 
radial harmonic trapping frequency 112 111 . The ID regime is 
realized when the transverse excitation energy hu± is much 
larger than the thermal energy ksT and chemical potential /i. 
While the YY thermodynamic equations do not directly yield 
the ID momentum distribution, n(k z ), here we show how they 
can be used to obtain the average kinetic energy per particle 



E kin /N = h 2 (k 2 )/2m, 



(1) 



where (k 2 ) 1 / 2 = [j 'dk z k 2 n(k z ) /N] 1 / 2 is the rms width of 
n{k z ), and N = J dk z n(k z ) is the total atom number. 

Solutions to the YY thermodynamic equations yield a 
unique value for the total energy per particle, E/N H, 
for a given temperature and interaction strength. Using 
the Helmann-Feynman theorem, the YY solutions can also 
be used to determine the local pair correlation function 



2 



m 0.4 




0.4- 



0.2 



- (b) 
















v//// 




1 o 5 y jfy 





10 10 10 10 10 10 10 

y 



Figure 1 . (Color online) (a) Solid lines show the kinetic energy per 
particle, E-^/N (in units of fcsT), from the exact YY solutions 
for a uniform ID Bose gas as a function of the chemical potential 
/i for different values of the dimensionless temperature parameter 
t. The gray dashed line is the ideal-Bose-gas result, corresponding 
to t — > oo (7 — ► 0). The crosses show the local kinetic energy per 
particle (without the contribution of the excited transverse modes) in 
the center of the trap for the experimental data presented in Fig. [3] 
(b) Same data as in (a), i.e., within the same window of values of 
(x/UbT, but presented as a function of the dimensionless interaction 
parameter 7 using the exact equation of state p = p(/i, T) @|, 



g^(0) Ella]- This gives the interaction energy per particle, 
Ei n t/N = \gpg^ 2 ' (0), and thus the kinetic energy per parti- 
cle is found from E kin /N = E/N - E int /N. 

The exact kinetic energy per particle for the uniform sys- 
tem, calculated for a wide range of interaction strengths and 
temperatures, is shown in Fig. Q] These results reveal that 
Ekin/N varies considerably as the degeneracy and interac- 
tion strength are changed; thus it is a useful observable for 
discriminating between different regimes of the ID Bose gas 
We identify three broad regimes of behavior. (;) For 
p < and \p\/ksT ^ 1, the system behaves as a classical 
(nondegenerate) ideal gas with E k i n /NksT — >• 1/2, as ex- 
pected from the classical equipartition theorem. As p — > 0~, 
the kinetic energy per particle, for sufficiently large values 
of t (negligible interactions), decreases as expected for a de- 
generate ideal ID Bose gas [gray dashed line in Fig. [TJ a )L 
E kin /Nk B T=lg 3/2 (\)/g 1/2 (\),whereg s (\) = j:Z 1 \ l /l s 
is a Bose function and A = e^/ kBT . (ii) Once p becomes pos- 
itive, the degenerate behavior is strongly affected by interac- 
tions. For 1 3> 1 (weak interactions, 7 <§; 1) the system is a qua- 
sicondensate and there is a significant reduction of E k i n /N 
from the classical equipartition value. (Hi) For t < 1 (strong 
interactions, 7 > 1) the system becomes fermionized, lead- 
ing to an increase in E k i n /N for both negative and positive p. 
For t < 1, and 1 < 7 <C t~ 1/2 fl] in the Tonks-Girardeau 
regime, E^/NksT becomes larger than the equipartition 



value of 1/2 |2 

In order to apply the YY results to trapped (nonuniform) 
quasi- ID systems as realized in experiment ! 1011 . we utilise the 
LDA. The trapping potential is U(r) = V (zj+'^muj^x 2 +y 2 ), 
where V{z) ~ \mui 2 z z 2 is the potential in the weakly confined 
longitudinal direction with lo z <C lj± i23ll . We calculate the 
particle number density p(z) by treating the trapped system 
as a collection of sufficiently small uniform systems of length 
Az, with the local chemical potential p(z) = p — V(z) ||8l 
lit] where p is the global chemical potential. The density of 
the trapped system in the ground transverse mode, po(z), is 
evaluated as po(z) — pyy[p(z), T], where Pyy[p, T] is the 
YY density for a uniform system. Similarly, the kinetic energy 
of the ground transverse mode is 



E^o = / dz£un[li.(z),T], 



(2) 



where £kin = E k - ln /Az is the kinetic energy density of a uni- 
form system of length Az. 

In the experiment lioll the system temperature was not suf- 
ficiently low that all transverse excitations were frozen out. 
We can account for this by observing that gpo -c hoj± so 
that the transverse excitations are well approximated as har- 
monic oscillator states with energies jfku± (j = 0, 1, 2, ...), 
where we have removed the zero-point energy ftid±. The 
transversely excited states can then be accurately described 
as independent ideal ID Bose gases with chemical potentials 
Pj{z) = p(z) — jhuj± ifioll . Accounting for the degeneracy 
factor j + 1, the ID position- and momentum-space densities 
for the atoms in transversely excited states are, respectively, 
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where (3 = l/k B T, and A T = (2irh 2 /niksT) 1 / 2 is the ther- 
mal de Broglie wavelength. Using Eq. © we obtain the ki- 
netic energy of the excited modes, £idn,e [cf. Eq. ([TJ]. Com- 
bining these results, the full YY description of the trapped 
quasi-lD system is given by the total density, p(z) = po(z) + 
p c (z), total atom number N = J p(z)dz, and the total average 
kinetic energy per particle E kin , t /N = (E kin . + E kin . c )/N. 

Before comparing the YY thermodynamic predictions for 
E k i n ,t/N with the experimental measurements of Ref. fioll . 
we outline the second theoretical approach used here — the 
SPGPE — which allows for the determination of both global 
thermodynamic quantities and the full momentum distribution 
in the weakly interacting regime. In the SPGPE approach, the 
system field operator for the lowest transverse mode is split 
into two parts, x &o(z) = ipc( z ) + i^iiz), representing coher- 
ent C0c) an d incoherent (tpi) regions lfl7tl . The coherent re- 
gion is defined by an energy cutoff e cut such that it contains 
all highly occupied modes. It can then be described as a clas- 
sical field (i.e. i/>c — > ipc) evolving according to the simple 
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Figure 3. (Color online) Temperature of a quasi- ID Bose gas in a 
harmonic trap (with cj z /2tt = 8.5 Hz and cj ± /2tt = 3280 Hz (H) 
as a function of the total kinetic energy iJkin.t determined by bi- 
modal Gaussian fits (circles) GH, SPGPE fits (triangles) |U, YY 
thermometry (crosses), and the classical ideal gas model (gray line). 
The inset shows the atom number (open circles, left axis) and -Ekin,t 
(right axis) as determined by bimodal Gaussian fits (open circles) 
llOn and SPGPE fits (triangles) to the momentum distributions as a 
function of the final rf evaporation frequency u T { , 



Figure 2. (Color online) Examples of the experimental momentum 
distributions (black dots), bimodal Gaussian fits (red dashed lines) 
and theoretical SPGPE best-fits (blue solid lines) yielding the tem- 
perature T and kinetic energy -Ekin.t- The bimodal Gaussian fits are 
difficult to distinguish from the SPGPE fits, but give slightly differ- 
ent values for the total kinetic energy (see inset of Fig. [3}. The ma- 
genta dot-dashed lines indicate the density of atoms in transversely 
excited modes n e (k z ). The finite experimental imaging resolution 
was modeled as a Gaussian with an rms width of 2.83 fim floll24ll . 
The uncertainty in the stated values of T and fi/ks is typically ±5 
nK (95% confidence interval). The values of v T i correspnd to final rf 
frequency for evaporative cooling. 



growth SPGPE 1171], 

# c =p{(^4)(/,-w 




(5) 



where C = — ■^d z + V(z)+g\^pc(z)\ 2 is the Gross-Pitaevskii 
operator, the parameters jd, M> ar, d T are the damping rate, 
chemical potential, and temperature of the reservoir, respec- 
tively, and dW is a complex Gaussian noise that is delta- 
correlated in time and space. This equation explicitly includes 
a projection (V) onto the coherent-region modes 12511 . and can 
be derived from a microscopic theory by tracing out the high 
energy modes that act as a thermal reservoir 11711 . In steady- 
state evolution the SPGPE samples ipc from a grand canonical 
density independent of the value of 7^. Thus the equilibrium 
density can be sampled in both position [pc( z ) = IV'C^)! 2 ] 
and momentum [nc(k z ) = \4>c{k z )\ 2 ] space, where the 
overline indicates time-averaging, and <f>c{k z ) is the spatial 
Fourier transform of ipc(z). 

The incoherent region (i.e., the longitudinal states of 
the ground transverse mode with low occupation) is well- 
described using the Hartree-Fock approximation 12611 . and has 



position density 
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dk z 
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+2gp c (z)-n] 



(6) 



where eg, = trk;/2m + V(z) |23fl. The incoherent re- 
gion momentum density n\(k z ) is obtained by a similar proce- 
dure, and the transverse ground mode distribution is no(k z ) = 
nc{k z ) + n\(k z ). Atoms in the excited transverse modes are 
treated as for the YY formalism, thus giving the full momen- 
tum distribution and total kinetic energy. 

Examples of best-fit momentum distributions obtained us- 
ing the SPGPE approach M27I1 are compared to the experimen- 
tal data l28ll in Fig. [2] We find quantitative agreement through- 
out the crossover from the nearly ideal Bose gas to the weakly 
interacting quasicondensate regime. From these fits we can 
determine both the kinetic energy per particle and the temper- 
ature for the data. 

In Ref. G3l the experimental momentum distributions were 
fitted with the sum of two Gaussians (see Fig. |2), giving the 
total number of atoms N, and the total kinetic energy -Ekin.t- 
The temperature T was obtained from the width of the broad- 
est Gaussian component using a classical ideal-gas model. We 
note that there is very little quantitative difference between the 
heuristic bimodal Gaussian momentum fits and those based on 
the microscopic SPGPE formalism. However, there is a dis- 
tinct difference in the temperatures extracted using the two 
methods (see Fig. [3]). 

We can also use the YY formalism for thermometry in this 
system, as there is a one-to-one correspondence between the 
temperature T and the kinetic energy per particle i?kin,t /N 
for a given atom number N. We compare the temperature es- 
timates of YY thermometry, the SPGPE momentum fits l29tl . 
and the broad Gaussian fits as a function of E^^/N in the 



4 



main panel of Fig. [3] A key result of this paper is that the YY 
temperatures collapse to the same curve as that obtained from 
the SPGPE fits to the full momentum distributions. 

The inset of Fig. [3] shows the total number of atoms as a 
function of the final rf frequency for evaporative cooling, as 
well as a comparison of -Ekin.t from the bimodal Gaussian fits 
and the SPGPE momentum fits. The temperature estimates 
from the SPGPE and Y Y methods differ for each final u> T { due 
to small differences between the SPGPE and bimodal Gaus- 
sian fits (see the comparison in Fig. |2). However, if the YY 
thermometry utilised the values of -Etin.t determined from 
the SPGPE fits, the temperature estimates would be identi- 
cal. Thus, the uncertainty in the temperature is due almost 
entirely to the experimental uncertainty in determining -Ekin,t- 
We note that the temperatures determined here t ypic ally lie 
above the estimates from the Gaussian fits of Ref. [1 1 Ofl — this 
illustrates the improved sensitivity of YY thermometry and 
SPGPE theory. 

The utility of YY thermometry is that it relies on a sin- 
gle measurement — the kinetic energy per particle — which can 
be readily obtained via focusing and straightforward density 
imaging. This provides a simpler alternative to thermome- 
try based on the measurement of density fluctuations of the 
gas J3| . It does not require any prior theoretical knowledge of 
the full momentum distribution, which is a challenging task in 
the strongly-correlated regime. 

Finally, the results we present in this paper provide further 
quantitative validation of the SPGPE in the regime of current 
experiments with ID quasi-condensates (see also H). As 
the SPGPE approach is suited to both equilibrium and dynam- 
ical simulations, this opens up an exciting avenue for explor- 
ing non-equilibrium phenomena in this system (e.g. quenches 
130D), which cannot be explored using the YY solutions or 
equilibrium quantum Monte Carlo techniques. 

In summary, we have studied the momentum properties of 
a finite-temperature ID Bose gas, and given exact results for 
the rms width of the momentum distribution. We have out- 
lined a procedure from which this simple quantity can be used 
for sensitive kinetic-energy thermometry of a quasi- ID har- 
monically trapped Bose gas using the YY formalism and the 
LDA. This method is applicable to all temperatures and re- 
pulsive interaction strengths. We have also performed a full 
characterization of the momentum distribution in the weakly 
interacting regime, and made a quantitative comparison with 
experimental data. Given the importance of momentum-space 
analysis in cold-atom research to date, our results provide a 
more complete picture of the homogeneous and trapped ID 
Bose gas, and open up new prospects for measuring and prob- 
ing these systems in the strongly-correlated regime. 
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